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Abstract. 

We introduce a Langevin equation describing the pinning-dcpinning phase 
transition experienced by Kardar-Parisi-Zhang interfaces in the presence of a bounding 
"lower- wall" . This provides a continuous description for this universality class, 
complementary to the different and already well documented one for the case of an 
"upper- wall" . The Langevin equation is written in terms of a field that is not an order- 
parameter, in contrast to standard approaches, and is studied both by employing a 
systematic mean-field approximation and by means of a recently introduced efficient 
integration scheme. Our findings are in good agreement with known results from 
microscopic models in this class, while the numerical precision is improved. This 
Langevin equation constitutes a sound starting point for further analytical calculations, 
beyond mean-field, needed to shed more light on this poorly understood universality 
class. 
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1. Introduction 

It was shown a few years ago that the introduction of a limiting or bounding wall 
into a Kardar-Parisi- Zhang (KPZ) interface model QjJ E] leads to quite different 
phenomenologies depending on whether the wall is an "upper" or a "lower" one 
[HI HI El Ej • This result, the origin of which can be traced back to the absence of height- 
inversion symmetry in KPZ interfacial dynamics Ej, has been verified for several 
discrete interfacial KPZ-like models [HI El El- In any of these cases, once a limiting 
wall is introduced, there are two different phases: a depinned one in which the KPZ- 
interface moves freely away from the wall, and a pinned one, with a finite stationary 
average distance from the wall. Separating these two phases there is a nonequilibrium 
phase transition, whose criticality encompasses also that of synchronizing extended 
systems 010], nonequilibrium wetting phenomena [TT H ITT j IT2]. transitions occurring in 
DNA alignment problems [SJ El, phenomena related to Burgers' turbulence [3|, bounded 
directed-polymers in random-media etc. Characterizing and distinguishing between 
the two above-mentioned pinning- depinning universality classes, with an upper or a 
lower wall respectively, is therefore a relevant task in many different contexts, as well 
clS cL chief theoretical problem. 

In terms of Langevin equations, an ideal framework to discuss universality issues, 
KPZ-like interfaces are described by the celebrated and profusely studied KPZ equation 



where h(r,t) is the interface height at position r and time t, A, D, and a are constants, 
and t] is a Gaussian white noise. In what follows, and without loss of generality, we 
take A = +D. Alternatively, we could also have fixed a given type of limiting wall, for 
instance, a lower, rigid substrate on top of which the interface grows, and observe the 
two different classes of depinning transitions under scrutiny depending on the sign of A 



Let us consider now equation (0) in the presence of an exponential, upper wall, i.e. 
including the additional term — 6exp(g h), with b > and q > 0. A transition takes 
place from a regime characterized by depinned interfaces, flowing to minus infinity 
for sufficiently negative values of a, to one with interfaces pinned to the wall, with 



HI: 



dh(r, t) 
dt 



a + A(V/i(r, t)) 2 + DV 2 h(r, t) + ar](r, t), 



(1) 
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exponentially cutoff positive values of h, above a certain threshold a c . This can also be 
visualized by performing a Cole-Hopf transformation, n = exp(h) (Tl|, which leads to, 

dn ^' - = DV 2 n(r, t) + an(r, t) - 6n(r, t) 1+<2 + tm(r, *Wr, t). (2) 
at 

This is a multiplicative noise equation (interpreted in the Stratonovich sense [TBI IT5]) 
defining a, well established by now, phase transition with a corresponding set of critical 
exponents that characterizes the so-called multiplicative noise 1 (MN1) universality class 
[H El El HZj • At the transition point and in the depinned phase the stationary average 
value of the order-parameter, n, is zero, while it is non-vanishing above the transition 
point. The critical exponents and other universal features in this class are not affected 
by the value of q, i.e. by the "impenetrability degree" of the wall; indeed, in microscopic 
models the wall is typically a hard substrate, corresponding to q —>■ oo. 

On the other hand, by introducing in equation ((H) a lower wall, bexp(-qh), 
hindering the interface height to take negative values, depinned interfaces flow towards 
plus infinity. In this natural order-parameter, equivalent to that for the upper- 

wall class, going to zero at the transition point, is m = exp(-h), and the corresponding 
depinning transition is in the so-called multiplicative noise 2 (MN2) universality class 

[HI El QUEUE]. 

Consequently, we perform the change of variables m = exp(-h) to obtain [T| 119j: 

dm [ T,t) = DV 2 m - 2D^¥- -am- 6m 1 " 9 + amr)(r, t), (3) 
at m 

where some space and time dependencies have been omitted to simplify the notation. 
Observe that owing to the (Vm) 2 /m term the equation becomes singular above the 
transition point, where fh — 0. Equation (jHJ) was studied in detail, both numerically 
and using mean- field approaches, in (7j, but no sound result could be obtained owing 
to the presence of the singular gradient term. Therefore, all the known results for 
this universality class come from numerical E] as well as some analytical (mean- 
field like) [THJ studies of discrete microscopic models. Let us also stress that direct 
numerical integrations of KPZ-like Langevin equations (before applying the Cole-Hopf 
transformation) are uncontrollable due to well-documented numerical instabilities [20 . 

Aimed at filling the gap between discrete and continuous levels of description for 
the MN2 universality class, in this article we show that the MN2 phenomenology can be 
captured by an alternative, multiplicative-noise, Langevin equation. To that purpose, 
we take the KPZ equation in the presence of a lower-wall and perform the change of 
variables, n = exp(h), leading to 

^ = DV 2 n + an + bn l ~ q + anr}(r, t), (4) 

which takes a particularly simple form for q — 1 although, as in the MN1 class, the 
precise value of q > is not expected to affect universal features. Note the remarkable 
difference between equations (j2j) and (J1J): while in the former n is an order-parameter 
for the MN1 transition, in equation (jlj) for the MN2 class, it is not, as it diverges at 
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the transition point and in the depinned phase. Therefore equation is not an order- 
parameter Langevin equation and it is m — 1/n that is to be monitored once the equation 
for n is integrated. Typically, Langevin equations are written in terms of a vanishing-at- 
the-transition order-parameter so that series expansions, truncation of power-series to 
lowest orders can be employed and the applicability of standard perturbative techniques 
is viable. This is in contrast with equation Furthermore, the noise amplitude 

diverges at the transition point. This apparent ill-behavior may be the reason why 
equation (J3J) has been ignored so far in the literature. 

In what follows we study the non-order parameter Langevin equation for the 
MN2 class, equation (jlj), using standard mean- field approaches and integrating it 
numerically. Despite of the presence of apparent pathologies and divergences at (and 
above) the transition point, we find that equation ([!} reproduces the previously known 
results for this universality class. In passing, we improve the numerical precision of 
the corresponding critical exponents. This constitutes a step forward in the general 
understanding of nonequilibrium phase transitions into absorbing states, allows for a 
better comparison with the MN1 class, gives a new starting point for future analytical 
approaches, and serves as an example of a phase transition best characterized in terms 
of an equation for a diverging field that is not an order-parameter. 

2. Mean-Field analysis 

Standard mean-field approaches usually neglect spatial and noise-induced fluctuations. 
For Langevin equations characterizing spatially extended systems with multiplicative 
noise, however, this approximation has proved to be too crude and unreliable in any 
dimension [Fj\ as both spatial structure and noise are relevant features. Therefore, more 
elaborated methods are required to obtain a sound qualitative picture of the transition 
[2*T| |22 EH]- A mean-field approach tailored to account properly for the noise term 
and the spatially- varying order-parameter consists in discretizing the Laplacian term as 
l/2dJ2j( n j ~ n i)i where d is the system dimensionality and the sum is over the nearest- 
neighbors of site i [221 E3j • When these latter are substituted by their average value 
n, a closed Fokker-Planck equation (which involves no approximation in the noise) is 
obtained for P(n, n). The stationary solution of such an equation is found by imposing 
the self-consistency requirement j2B] 



Note that this approach preserves the crucial role played by the multiplicative noise and 
includes the spatial coupling even if in an approximate (self-consistent) way. A detailed 
discussion of the results obtained with this method for the MN1 class can be found in 
[21]. Applying this procedure to equation (JU) one readily obtains 



n = 



Jq 00 dn n P(n, n) 



(5) 



Jq 00 dn P(n, n) 




(6) 
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and after defining I p (n) = / dn n p P(n,n), 
- h{n) _ /-i(ra) 

n= TFv m = r / — \ * (7) 
J (n) J (n) 

To evaluate the scaling behavior of I p near the transition, when n becomes large, we 
first do the substitution z = 2Dn/a 2 n and then expand the newly generated term 
exp{— 2b/(a 2 ( 1 ~ qS> q) x (z/2Dn) q } to first order to obtain 

I p (n) « A p n?+i + B v fi^~\ (8) 

where 7 = 2 (a — D)/a 2 and Ap, 5 p are expressed using the Gamma function as 

{2D\ y+p . 

A p={jr) r(- P - 7 ), 

/2m p+7 - g 2b , . 

B r =-(-^) ^r(-p- 7 + 5) . ( 9) 

A direct calculation then yields n ~ |a + cr 2 /2\ and m ~ \a + a 2 /2\ l / q . Therefore, 
the order-parameter critical exponent is (3 = 1/q [23]. 

It is instructive to compare these results with those obtained using the same 
technique for equation ©, (3 = max (1/q, a 2 /{2D)) j22H2H- That is, in the MN1 case 
two possible values of (3, reminiscent of the strong- and weak-coupling regimes of the 
KPZ dynamics [2J, appear already at this self-consistent mean-field level. In contrast, for 
equation (@J), there is no strong-coupling regime, which would be characterized by a non- 
universal, noise- amplitude-dependent, value of (3. This fact is related to the presence 
of a single cut-off in the stationary probability distribution, equation while two 
different cut-offs, and therefore two different mechanisms controlling the scaling, appear 
for equation (j2J) |2I]- The implications of this property, as well as its connections with 
the high-dimensional behavior of the KPZ dynamics, will be analyzed elsewhere. 



3. Numerical integration of stochastic differential equations 

In order to integrate equation (j3J) as efficiently as possible we have employed a recently 
introduced split-step scheme for the integration of Langevin equations with non-additive 
noise. In this scheme, the Langevin equation under consideration is studied on a 
lattice and separated in two parts: the first one includes only deterministic terms and 
is integrated at each time-step using any standard integration scheme: Euler, Runge 
Kutta, etc [23] (here we have chosen a simple Euler algorithm). The output of this step 
is used as the input to integrate (along the same time-step) the second part which may 
include the linear deterministic term and, more importantly, the noise. This is done by 
sampling in an exact way the probability distribution function associated with this part 
of the equation. In the case under study (noise proportional to the field) the second 
step corresponds to sampling a log-normal distribution solution of d t n = an + anrj (for 
more details see |2H! and |2Z1)- Therefore the two-step algorithm is finally specified by: 

m(i, t) = n(i, t) + (b n(i, t) 1_<? + D[n(i + l,t) + n(i - 1, t) - 2n(i, t)]) dt (10) 
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and, 



n(i, t + dt) — rii(i, t) exp (adt + a Vdt r\ 



where i] is a random variable extracted from a Normal distribution with zero-mean and 
unit variance. Note that the linear deterministic term can be included in either the first 
or the second step, or partially incorporated in both of them. 

We have considered one-dimensional lattices, and fixed a = 1, b = 1, D = 0.1, 
space-mesh dx = 1, and time-mesh dt = 0.1 (note that in this scheme dt can be 
taken larger than in usual integration algorithms |27|). As initial condition we take 
n(r, t = 0) = 3. We take q = 4 for all simulations except for results presented in figure 
a) where we show that asymptotic results do not depend on the value of q, as long 
as q > 0. We then iterate the dynamics by employing the previous two-step integration 
algorithm, using parallel updating, at each site i. 

A summary of our main findings follows. First, to accurately determine the critical 
point, we perform decay experiments and average over many independent runs in a 
system of size L = 2 17 . At criticality, a c = —0.143668(3), the average density, rh = (1/n) 
decays as a power-law with an associated exponent 9 ~ 0.229(5) (see FigJIJ. This is 
to be compared with the previous estimates 9 ~ 0.215(15) [7] and 9 ~ 0.228(5) [12J. 
On the other hand, for smaller system-sizes, we observe saturation at this value of a c , 
and the scaling of the saturation values for different system sizes (inset Fig^a)) gives 
P/v w 0.335(5) (to be compared with 0.34(2) in 0). 




Figure 1. A: Log-log plot of the time decay of the average order-parameter at the 
critical point, a c = -0.143668, for system sizes 2 6 , 2 7 , 2 8 , 2 9 , 2 10 , 2 U , and 2 17 
respectively and q = 3. In the inset, the average saturation values of the previous 
curves are plotted as a function of the system-size, L, in double-logarithmic scale. B: 
Averaged order-parameter decay for equation (@J with q — 0.5, 1, 2, 4 (from bottom to 
top) in a system of size L = 2 17 . For any q > we observe the same asymptotic decay 
exponent at the same critical point location. 



For other values of q we have verified that, as shown in figure 2, none of the 
previously reported exponents, nor the location of the critical point, are altered, 
although q— dependent transient effects exist. The invariance of a c against changes 
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in q can be understood from the fact that a c corresponds to the value of a for which 
depinned interfaces, arbitrarily far from the wall, invert their direction of motion; this 
is not affected by the nature of the wall, i.e. by the value of q. Additionally, the order- 
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Figure 2. A: Log-log plot of the saturation values of m for different values of a 
nearby the critical point. From the slope we estimate (3 — 0.332(5). B: Log-log plot of 
the time growth of the averaged height h = — log(n) for the same value of a and same 
system sizes as in figure 1. In the inset, the average saturation values of the previous 
curves are plotted as a function of the system-size, L, in a log-log scale. 



parameter exponent (3 ~ 0.332(5) has been measured for L = 2 17 (see also figure 2). 
The previous best estimation was j3 = 0.32(2) jjj. 

We have also studied the scaling properties of the height field h = — log(n) (see 
Fig|2fb)). An analysis analogous to that presented above for m leads to h{t) m i°- 323 ( 10 ) ; 
h ~ \a — a c |~ a48 ( 3 ), and h(a c ,L) ~ L 0A8( - 2 \ which define a set of critical exponents for 
h analogous to those for m : Q h w -0.323(10), (3 h w -0.48(3), and (3 h /v w -0.48(2) 
(see figure 2b). The previous best estimations for these exponents were: —0.355(15), 
—0.52(2) and, —0.52(2) respectively j7j. The minus signs just reflect the fact that h 
diverges at the transition point. A summary including our best estimates for the critical 
exponents can be found in table 1. Finally, we stress that our results are compatible 
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P 


P/v 


V 


z 


m 


0.229(5) 


0.332(5) 


0.335(5) 


0.99(3) 


1.46(5) 


h 


-0.323(10) 


-0.48(3) 


-0.48(2) 


1.0(1) 


1.48(10) 



Table 1. Critical exponents for the MN2 class in d = 1 as calculated in this paper. 
Scaling of saturation times gives z which is consistent with the value z = (3/{9v) 



with the theoretical predictions z = 3/2 and v — 1, expected to be valid for both MN1 
and MN2 jUO]- Also, all the standard scaling laws between exponents are satisfied j^j. 
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4. Discussion 

We have studied the dynamics of KPZ-like interfaces with A > bounded by a lower- 
wall. The results and phenomenology differ from that of upper-walls. By performing 
a Cole-Hopf or logarithmic transformation, the resulting order-parameter Langevin 
equation (JHJ) is singular and no satisfying result can be derived from it. Instead, the 
main result of this paper, is that a sound Langevin equation (J3J) can be written in terms 
of a non-order-parameter field which diverges at the transition. For such an equation we 
have performed (i) a self-consistent mean-field analysis leading to the result (3 = 1/q, and 
no trace of any strong-coupling regime (noise-dependent j3 exponent value) contrarily 
to what happens for the upper-wall case; (ii) extensive numerical integrations of the 
stochastic equation using a recently introduced very-efficient numerical scheme. The 
obtained critical exponent values are in good agreement with previously known ones 
measured in simulations of discrete models, and improve the level of accuracy and 
precision. 

In summary, we have shown that an apparently ill-behaved non-order-parameter 
Langevin equation constitutes a sound continuous representation of the pinning- 
depinning transition experienced by interfaces in the Kardar-Parisi- Zhang class under 
the presence of a bounding lower-wall. Performing further analytical, renormalization 
group analyses of the present Langevin equation remains as a challenging task. 
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